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Abstract 

This  paper  theoretically  studies  the  water  transport  phenomena  in  PEM  fuel  cells,  mainly  investigating  the  transient  behavior  in  the  gas  diffusion 
layer  (GDL),  catalyst  layer  (CL)  and  proton  exchange  membrane  (PEM).  In  the  PEM,  both  diffusion  and  electro-osmosis  processes  are  considered, 
while  in  the  GDL  and  CL,  only  diffusion  process  is  taken  into  account.  The  process  of  water  uptake  is  employed  to  account  for  the  water  transport 
at  the  interface  between  the  PEM  and  CL.  The  results  indicate  that  the  water  content  in  the  PEM  and  the  time  for  reaching  the  steady  state  in  the 
start-up  process  are  influenced  by  the  humidification  constant,  k,  the  humidification,  and  the  thickness  of  PEM.  The  rise  of  the  k  increases  the  water 
content  in  the  membrane  and  shortens  the  time  for  reaching  the  steady  state.  Insufficient  humidification  causes  relatively  small  water  content  and 
long  steady  time.  When  the  PEM  is  thinner,  the  water  is  more  uniformly  distributed,  the  water  content  gets  higher,  and  the  time  for  reaching  the 
steady  state  is  distinctly  shorter. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Water  management  gets  to  be  a  crucial  problem  in  the  recent 
research  in  PEM  fuel  cells,  in  which  adjusting  associated  param¬ 
eters  is  an  important  consideration.  Okada  demonstrated  a  series 
of  studies  on  water  transport  and  distribution  in  the  electrolyte 
membrane  of  fuel  cells  to  investigate  the  influence  of  the  mem¬ 
brane  characteristics,  operating  conditions,  and  contamination. 
An  one-dimensional  mathematical  model  involving  water  dif¬ 
fusion  and  electro-osmotic  process  in  a  membrane  to  discuss 
the  water  content  profile  under  various  operating  conditions  was 
proposed  by  Okada  et  al.  [1].  The  results  indicated  that  the  water 
content  profile  was  greatly  influenced  by  the  electrical  current 
density  and  water  transfer  coefficient.  In  1998,  Okada  et  al.  [2,3] 
experimentally  and  analytically  investigated  the  effects  of  con¬ 
tamination  on  water  transport  in  a  half  membrane  of  the  anode 
side.  They  found  that,  when  the  PEM  was  polluted,  the  water 
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molecules  were  dragged  by  the  contaminator  rather  than  hydro¬ 
gen  ions,  which  caused  the  drop  of  electrical  current  density. 
Okada  et  al.  [4]  extended  the  previous  model  to  the  whole  mem¬ 
brane,  and  conducted  some  calculations  according  to  the  actual 
experimental  conditions.  Later,  Okada  et  al.  [5,6]  theoretically 
analyzed  the  influence  of  contamination  on  the  whole  mem¬ 
brane  characteristics  in  the  conditions  of  contaminated  anode 
input,  and  the  cathode  input  with  finite  thickness  of  contami¬ 
nated  region. 

Eikerling  et  al.  [7]  presented  a  model  to  predict  the  water 
transport  process  under  the  electro-osmosis  drag  and  capil¬ 
lary  force.  Bernardi  and  Verbrugge  [8]  put  forward  a  model  to 
describe  the  physical  process  in  the  PEM,  CL  and  GDL  on  the 
cathode  side,  mainly  investigating  the  cathode  polarization,  the 
water  transport  in  the  PEM  and  CL,  and  the  utilization  ratio  of 
the  CL.  Bernardi  and  Verbrugge  [9]  examined  the  mass  transport 
and  reaction  process  in  the  region  including  the  anode  CL,  PEM, 
cathode  CL  and  cathode  GDL.  Singh  et  al.  [10]  further  extended 
the  one-dimensional  model  of  Ref.  [9]  into  two-dimensional 
model  to  analyze  the  dimensional  effects.  It  was  pointed  out  that 
the  voltage  resulted  from  the  two-dimensional  model  was  appre- 
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Nomenclature 

a  water  activity 

C[j  water  concentration  in  each  layer  (mol  cm-3) 

Dm  water  diffusivity  in  PEM  (cm2  s-1) 

D(0)  the  constant  term  in  Taylor  series  of  water  diffu¬ 
sivity  in  PEM  (cm2  s-1) 

D(1)  the  first  order  term  in  Taylor  series  of  water  dif¬ 
fusivity  in  PEM  (cm2  s-1) 
water  vapor  diffusivity  in  gas  diffusion  layer 
(cm2  s-1) 

D ^  water  vapor  diffusivity  in  catalyst  layer  (cm2  s- 1 ) 

F  Faraday’s  constant  (C  mol-1) 

I  electrical  current  density  (A  cm-2) 

7h2o  water  flux  (mol  cm-2  s  - 1 ) 
k  humidification  constant  (cm  s  - 1 ) 

Mw  molecular  weight  of  water  (g  mol-1) 

P  pressure 

f  mean  pore  size  (cm) 

R  universal  gas  constant  ( J  mol  K- 1 ) 

w  water  transfer  coefficient 

u/01  the  constant  term  in  Taylor  series  of  water  transfer 
coefficient  in  PEM 

the  first-order  term  in  Taylor  series  of  water  trans¬ 
fer  coefficient  in  PEM 

Greek  letters 

sc  porosity  of  catalyst  layer 

£d  porosity  of  gas  diffusion  layer 

£  electrochemical  reaction  rate  (A  cm-3) 

fj  molar  fraction  of  specie  j 

X  water  content  in  membrane  (mol  H20/mol  SO3  -) 

Subscripts 

c  catalyst  layer 

d  gas  diffusion  layer 

g  gas 

1  liquid 

m  proton  exchange  membrane 


ciably  less  than  that  from  the  one-dimensional  model.  Springer  et 
al.  [  1 1  ]  proposed  a  one-dimensional  isothermal  model  to  analyze 
the  water  transport,  and  found  that  the  ratio  of  water  molecule 
flux  to  hydrogen  ion  flux  was  quite  less  than  the  directly  mea¬ 
sured  electro-osmosis  coefficient.  Springer  et  al.  [12]  further 
indicated  that,  when  the  cathode  fuel  was  air,  the  potential  loss 
of  PEM  was  only  a  small  part  of  the  overall  loss,  which  mainly 
occurred  in  the  cathode.  Fuller  and  Newman  [13]  theoretically 
investigated  the  water  and  thermal  management  in  fuel  cells. 
Bradean  et  al.  [14]  presented  a  model  of  a  half  fuel  cell  on  the 
cathode  side  to  discuss  the  influence  of  operating  conditions  on 
the  cell  performance.  They  indicated  that  the  diffusion  process 
was  the  main  mode  of  the  mass  transport,  and  the  occurrence  of 
super  saturation  of  water  vapor  was  definitely  associated  with 
the  cathode  thickness  and  the  reaction  rate. 


Nguyen  and  White  [15]  developed  a  one-dimensional  model 
to  investigate  the  water  transport  in  the  membrane  under  the 
diffusion  and  electro-osmosis  processes,  coupled  with  the  heat 
transfer  between  the  solid  and  gas  phases,  involving  the  latent 
heat  absorbed/emitted  in  evaporation/condensation.  The  results 
showed  that,  when  the  current  density  was  over  1  A  cm-2,  the 
Ohmic  loss  in  the  membrane  played  an  important  role  in  the 
overall  potential  loss.  Yi  and  Nguyen  [16]  improved  the  model 
in  Ref.  [15]  to  involve  the  convection  process  induced  by  the 
pressure  gradient,  and  analyzed  the  temperature  profile  along  the 
flow  passages.  Later,  He  et  al.  [17]  proposed  a  two-dimensional 
model  to  study  the  two-phase  flow  phenomena  in  fuel  cells.  To 
clarify  the  water  effects,  Yan  et  al.  [18]  examined  the  effects  of 
flow  distributor  geometry  and  diffusion  layer  porosity  on  reac¬ 
tant  gas  transport  and  performance  of  PEM  fuel  cells.  In  Ref. 
[18],  the  porosity  of  gas  diffusion  layer  was  varied  to  account 
the  effects  of  water  content.  They  found  that  the  porosity  of  gas 
diffusion  layer  has  significant  impact  on  the  cell  performance. 

The  water  content  in  the  membrane  affects  the  electrical  resis¬ 
tance  to  influence  the  cell  performance.  An  experimental  study 
on  Nation  117  membrane  performance  on  water  transport  at 
30 °C  was  conducted  by  Zawodzinski  et  al.  [19].  Hinatsu  et 
al.  [20]  also  performed  an  experiment  to  investigate  the  water 
uptake  process,  and  obtained  the  performance  of  membrane  sub¬ 
merged  in  the  liquid  water  at  a  temperature  from  25  to  130  °C. 
Futerko  and  Hsing  [21]  deduced  a  model  to  predict  the  water 
content  variation  in  the  membrane,  and  compared  the  results 
with  the  available  experimental  results.  Recently,  the  dynamic 
performance  of  PEM  fuel  cell  was  examined  by  Chen  et  al. 
[22].  In  their  work,  the  transient  behaviors  of  the  water  transport 
across  the  membrane  were  examined,  as  well  as  the  influences  of 
several  physical  parameters  on  the  characteristic  time  to  reach 
the  steady  state.  The  detailed  thermal  and  water  management 
in  the  membrane  of  PEM  fuel  cells  is  investigated  numerically 
by  Yan  et  al.  [23].  The  coupling  effects  of  mass  diffusion  and 
temperature  gradient  on  the  water  distribution  of  the  membrane 
were  taken  into  account  with  consideration  of  the  temperature- 
dependent  diffusivity. 

The  literature,  as  described  above,  was  mainly  restricted  in 
the  steady- state  performance  and  mass  distribution  in  the  fuel 
cell,  but  lacks  the  transient  analysis.  The  present  work  aims  to 
investigate  the  transient  process  to  obtain  the  time  for  reaching 
the  steady  situation  and  to  understand  the  evolution  of  water 
content  under  various  conditions.  A  one-dimensional  transient 
model  is  further  developed  by  adding  the  mass  transport  pro¬ 
cess  in  GDL  and  CL  into  the  steady  model  derived  by  Okada’s 
research  group,  which  only  considered  the  transport  process  in 
the  membrane. 

2.  Theoretical  analysis 

2.1.  Assumptions 

In  the  present  work,  a  systematically  parametric  study  is  con¬ 
ducted  to  examine  the  transient  behavior  of  the  water  content 
across  the  GDL,  CL  and  membrane  under  the  influences  of  var¬ 
ious  physical  parameters.  A  special  attention  is  paid  to  the  time 
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needed  for  the  water  distribution  across  the  system  to  reach  the 
steady  state,  denoted  as  ts,  under  influence  of  different  physical 
parameters.  Since  the  thicknesses  of  the  PEM,  GDL  and  CL  are 
quite  small  comparing  to  the  length  and  width,  the  fuel  cell  is 
regarded  as  a  one-dimensional  system.  For  convenience,  other 
assumptions  are  made  as  the  following: 

(1)  The  system  is  one-dimensional  and  isothermal. 

(2)  In  the  CL  and  GDL,  the  diffusion  is  the  only  mass  transport 
mode. 

(3)  The  membrane,  CL  and  GDL  are  homogeneous. 

(4)  The  electrochemical  reaction  only  occurs  in  the  CL  on  the 
cathode  side. 

(5)  In  the  GDL  and  CL,  water  only  exists  in  vapor  phase,  and 
could  be  supersaturated. 


mulations  and  relevant  physical  parameters  in  this  section  for  the 
convenience  of  following  discussions. 

As  the  water  transport  in  the  PEM  is  induced  by  the  diffusion 
and  electro-osmosis  processes,  the  water  flux  in  the  PEM  is 
described  as  the  following: 

t  r»m  ^C3,l  |  I  ^ 

AbO  =  —D  — - b  —w  (1) 

ox  F 

where  Dm,  c33,  /,  and  w  denote,  respectively,  the  water  diffusiv- 
ity  coefficient,  water  content,  current  density,  and  water  transfer 
coefficient  of  the  PEM  due  to  the  electro-osmosis  drag.  And  F 
is  Faraday’s  constant.  The  water  content  evolution  in  the  mem¬ 
brane  is  expressed  as  below: 


dc^l  dJu20  d_  ( jynfc 33  _ 
d  t  dx  dx  \  dx  F  ) 


(2) 


2.2.  Governing  equations 

The  physical  model  is  schematically  shown  in  Fig.  1.  The 
origin  is  set  at  the  border  of  the  anode  GDL,  which  is  away  from 
the  interface  between  the  anode  GDL  and  CL  at  a  distance  of  A, 
away  from  the  interface  between  the  anode  CL  and  PEM  at  B, 
away  from  the  interface  between  the  cathode  CL  and  PME  at  C, 
away  from  the  interface  between  the  cathode  CL  and  GDL  at  D, 
and  away  from  the  border  of  cathode  GDL  at  E.  The  mass  trans¬ 
port  equations  in  the  GDL,  CL  and  PEM  are  expressed  below. 
Although  the  theoretical  model  is  essentially  an  extension  from 
Okada  [1-7],  we  still  present  the  details  of  the  mathematical  for- 


Since  the  water  diffusivity  and  water  transfer  coefficient  vary 
with  water  content  in  the  membrane,  they  are  expressed  in  the 


forms  of  Taylor  series: 

Dm  =  D<0)  +  D(1)C3,!  +  D(2)c|  ,  H - 

(3) 

W  =  u/0)  +  +  ub-’cj  ]  +  •  ■  ■ 

(4) 

Substituting  Eqs.  (3)  and  (4)  into  Eq.  (2)  yields: 


dc3,l 

dt 


D(0)  d 

a2* 


—  w 
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dx 
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d2C3,l 

d2x 


-hD 


(1) 


f  dc3,i X  2 

V  dx 
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Cathode 


21  (2)  dc3  J 

—  WK  ’c3  1 - 

F  ’  dx 


+ 


(5) 


In  Eq.  (5),  the  order  of  terms  in  the  first  pair  of  square  bracket 
is  much  higher  than  that  in  the  second  one,  thus  the  terms  in 
the  second  square  brackets  are  neglected.  Then  the  equation  is 
reduced  into  the  form: 


(6) 


According  to  Eq.  (6),  the  water  diffusivity  can  be  assumed 
constant,  and  the  water  transfer  coefficient  should  include  the 
constant  and  first  order  terms  of  the  Taylor’s  series. 

In  the  GDL,  only  diffusion  process  is  taken  into  account  to 
describe  the  water  vapor  transport.  As  the  mean  pore  size  in  the 
GDL  is  far  larger  than  the  mean  free  path  of  the  water  vapor 
molecules,  it  is  appropriate  to  express  the  transport  process  as  a 
binary  diffusion  process.  Then  the  water  vapor  flux  in  the  GDL 
is  described  as  below: 


Jw2o  =  -^wVG,  g 


(7) 


0 


AB 


CD 


where  D\  is  calculated  from  the  following  equation: 


GDL 


PEM  GDL 


Fig.  1 .  Schematic  diagram  of  the  physical  system. 


(8) 
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where  Sd  is  the  porosity  of  the  GDL,  fj  the  molar  fraction  of 
specie  j,  and  ■  is  the  binary  diffusion  coefficient  of  water 
vapor  and  specie  j. 

Substituting  Eq.  (8)  into  the  mass  conservation  equation 
yields  the  diffusion  equations  of  water  vapor  in  the  GDL: 


£d 


9cpg 

dt 


8 

dx 


nb  dcbg 
w  dx 


Sd 


^c5,g 

dt 


d 

dx 


^c5,g 

dx 


0  <  x  <  A 

D  <  x  <  E 


(9) 

(10) 


where  c i,g  and  C55g  are  the  water  vapor  content  in  the  anode  GDL 
and  cathode  GDL,  respectively. 

Similarly,  only  diffusion  process  is  considered  in  describing 
the  water  vapor  transport  in  the  CL.  As  the  pore  size  in  CL  is 
relatively  small,  the  collision  mainly  occurs  between  the  water 
molecules  and  the  wall,  which  is  quite  different  from  the  sit¬ 
uation  in  GDL.  Therefore,  Knudsen  diffusion  is  employed  to 
describe  the  diffusion,  and  then  the  water  vapor  flux  in  the  CL 
is  obtained: 


Ju2  o  =  —  g  (11) 

where  Knudsen  diffusion  coefficient  can  be  calculated  by  the 
following  expression: 


2.4.  Boundary  conditions 


The  relative  humidity  is  specified  as  constant  in  the  two  outer 
borders  of  the  physical  model.  In  the  interfaces  between  the  GDL 
and  CL,  the  water  concentration  and  water  flux  are  continuous. 
In  the  interface  between  the  CL  and  the  PEM,  the  water  transport 
is  treated  with  the  water  uptake  process.  The  water  content  in  the 
PEM  is  relevant  to  the  water  activity  of  the  surrounding  water 
vapor,  which  is  defined  as  below: 


a  = 


(16) 


where  c|at  represents  the  saturated  water  vapor  concentration 
under  a  certain  temperature.  The  water  content  in  the  membrane 
is  obtained  by  sorption  isotherm  A  (a): 

X  =  —  =  A{a)  (17) 

cH+ 

where  cH+  is  the  proton  molar  concentration  in  the  membrane, 
and  A(a)  is  calculated  by  the  following  empirical  formula,  which 
is  derived  from  the  experimental  results  for  Nafion  1 17  at  80  °C: 


A  =  0.3  +  6a[l  —  tanh  (a  —  0.5)]  +  3.9  \[a 


1  +  tanh 


a-  0.89V 
0.23  J  _ 


(18) 


(12) 


Generally,  the  boundary  conditions  of  the  model  are  listed  as 
below: 


where  x/r  is  the  correction  factor,  f  the  mean  pore  size,  Mw  the 
water  vapor  molecular  weight,  and  R  is  the  universal  gas  con¬ 
stant. 

Letting  the  mass  conservation  replaced  by  Eq.  (11),  it 
becomes: 


Sc 


dc2,g 

dt 


d 

dx 


A  <  x  <  B 


(13) 


x  —  0  :  Q,g  —  Qn 
x  =  A  :  ci,g  =  C2,g 


x  =  A  : 


b  dcPg  dc2,g 

w  dx  w  dx 


x  =  B  : 


dcfg 

dx 


=  k[Ca 


C3,l\ 


(19) 

(20) 

(21) 

(22) 


It  is  presumed  that  the  water  vapor  generated  from  the  electro¬ 
chemical  reaction  is  uniformly  distributed  in  the  cathode  CL. 
Then  the  diffusion  equation  becomes: 


Sc 


d^4,g 

dt 


d 

dx 


^c4,g 

dx 


C  <  x  <  D 


(14) 


where  sc  is  the  porosity  of  the  CL,  C2,g  and  c^g,  respectively, 
denote  the  water  vapor  concentrations  in  the  anode  CL  and  cath¬ 
ode  CL,  and  £  is  the  electrochemical  reaction  rate. 


2.3.  Initial  conditions 


x  =  B  :  k[Ca  -  c3,i]  =  -D(0)^  +  f  [«/0)  +  ^(1)c3,i] 

ox  F 


(23) 


x  =  C:  k[Cc  —  c3  i]  =  D(0)  ^  f  [u/o>  +  wWC3  ] 

ox  F 


(24) 


x  =  C  :  =k[Cc-cxx 


w  dx 
x  =  D  :  c4,g  =  c5, 


(25) 

(26) 


The  initial  water  concentration  profiles  in  the  PEM,  CL  and 
GDL  are  assumed  to  be  the  steady-state  distributions  in  the  pre¬ 
vious  operation,  cf  •(*): 


x  —  D  \ 
x  =  E  : 


k  ^4,g  _  b  dc5,g 

w  dx  w  dx 

c5,g  —  Qn 


(27) 

(28) 


t  =  0  :  Cij(x,  0)  =  cf(x)  (15) 

when  i  =  1 , 2, 4, 5 ,j  =  g,  corresponding  to  the  state  in  the  GDL  and 
CL.  When  i  =  3 ,j  =  1,  corresponding  to  the  state  in  the  membrane. 


where  cm  denotes  the  specified  water  concentration,  k  the  rate 
constant  of  Henry’s  law  accounting  for  the  water  flux  entering  or 
exiting  the  membrane  driven  by  the  water  concentration  gradient 
in  the  PEM,  which  is  also  called  the  humidification  parameter. 


Table  1 

Values  of  parameters  used  in  this  work 
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Parameter 

Symbol 

Value 

Pressure  at  anode  side  (atm) 

P  a 

3 

Pressure  at  cathode  side  (atm) 

Pc 

3 

Cell  temperature  (°C) 

Tcell 

80 

Thickness  of  GDL  (cm) 

Ld 

0.03 

Thickness  of  CL  (cm) 

Lc 

0.002 

Thickness  of  PEM  (cm) 

Em 

0.01 

Effective  porosity  of  GDL 

£d 

0.4 

Effective  porosity  of  CL 

£C 

0.3 

Effective  porosity  of  PEM 

£m 

0.28 

The  constant  term  in  Taylor  series  of  water  diffusivity  in  PEM  (cm2  s~ 1 ) 

T>(°) 

1.25  x  10-5 

The  constant  term  in  Taylor  series  of  water  transfer  coefficient  in  PEM 

w(0) 

0.0 

The  first-order  term  in  Taylor  series  of  water  transfer  coefficient  in  PEM 

u/1) 

128.32 

Water  vapor  diffusivity  in  gas  diffusion  layer  (cm2  s“ 1 ) 

2.13  x  10“5 

Water  vapor  diffusivity  in  catalyst  layer  (cm2  s~ 1 ) 

1.39  x  10-5 

Ca  and  Cc,  respectively,  represent  the  liquid  water  concentra¬ 
tions  at  the  interfaces  between  the  PEM  and  anode  CL,  and 
between  the  PEM  and  the  cathode  CL,  which  are  at  equilibrium 
with  the  liquid  water  concentrations  in  the  CL.  Ca  and  Cc  are 
calculated  from  the  following  expressions: 


Cu  —  Tp[+ 


c2,g 


x=B 


■sat 

w 


(29) 


Cc  —  d"p[+  .A 


c4,g  x=c 


■sat 

'W 


(30) 


3.  Results  and  discussions 


Based  on  the  theoretical  model  with  boundary  and  initial 
conditions  shown  in  the  above  section,  we  perform  a  system¬ 
atically  parametric  study  to  examine  the  transient  behavior  of 
the  water  content  across  the  GDL,  CL  and  membrane  under 
influences  of  various  physical  parameters.  These  parameters 


include  the  current  density,  i,  the  initial  current  density,  io, 
the  relative  humidity  of  water  vapor  on  the  cathode  or  anode 
side,  the  thickness  of  PEM,  and  the  humidification  parameter, 
k.  The  values  of  these  parameters  employed  in  the  theoreti¬ 
cal  model  are  corresponding  to  the  base  case  at  Tcen  =  80oC, 
P-  3  atm,  Lm  =  0.01  cm,  io  =  0  A  cm-2,  &  =  0.01  cm  s-1,  and  rel¬ 
ative  humidity  100%.  When  the  physical  effect  of  certain  param¬ 
eter  is  evaluated,  the  parameter  varies  with  the  others  fixed  the 
same  with  the  base  case.  The  parameters  used  in  this  work  are 
shown  in  Table  1 . 

Lig.  2  presents  the  steady-state  distributions  of  the  water  con¬ 
tent  under  four  different  operating  current  densities,  with  the 
humidification  parameter  k  0.01  cm  s_1 ,  the  relative  humidity  of 
water  vapor  on  the  cathode  and  anode  sides  100%,  and  the  thick¬ 
ness  of  PEM  0.01  cm.  When  a  content  difference  exists  between 
the  PEM  and  its  outside,  k  denotes  the  magnitude  of  the  equi¬ 
librium  rate.  The  left  figure  in  Fig.  2  shows  the  distributions  of 
water  vapor  content  in  the  anode  GDL  and  CL,  the  middle  one 
shows  the  distribution  of  liquid  water  content  in  the  PEM,  and 


( x-L  i -L  2~L3 )/( L4+L5 ) 


Fig.  2.  The  steady  distributions  of  the  water  content  under  four  different  operating  current  densities,  where  the  humidification  constant  k  is  0.01  cm  s  1 ,  the  relative 
humidity  of  water  vapor  on  the  cathode  and  anode  sides  100%,  and  the  thickness  of  PEM  100  pan. 
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the  right  one  indicates  the  distributions  of  water  vapor  content  in 
the  cathode  GDL  and  CL.  Inside  the  PEM,  the  water  content  near 
the  anode  side  is  relatively  low  due  to  the  electro-osmotic  effects. 
As  the  current  density  i  increasing,  the  electro-osmotic  effects 
are  enhanced,  and  the  water  content  difference  between  the  two 
sides  of  the  PEM  becomes  more  evident.  Since  the  water  content 
near  the  anode  side  of  PEM  is  relatively  low,  the  water  vapor 
content  in  anode  GDL  and  CL  will  be  reduced  to  compensate 
the  water  depletion  in  PEM.  The  phenomena  in  the  cathode  are 
just  opposite.  The  water  content  near  the  cathode  side  of  PEM 
is  relatively  high,  and  the  electrochemical  reaction  in  cathode 
CL  generates  liquid  water.  Therefore,  the  water  content  in  cath¬ 
ode  CL  is  higher  than  that  of  the  cathode  GDL.  In  a  separate 
numerical  simulation,  it  is  found  that  when  the  humidification 
parameter  k  is  small,  the  flux  of  the  water  vapor  absorbed  by  the 
PEM  is  low,  resulting  high  the  water  content  in  anode  CL. 

Fig.  3(a)  and  (b)  presents  the  evolution  of  the  water  content  in 
PEM  from  steady  distribution  at  z'o  =  0.1  A  cm-2  with  humidifi¬ 
cation  constant  being  0.01  and  0.001  cm  s_1 ,  respectively,  when 
the  operating  current  jumped  to  i  =  0.1  A  cm-2.  It  is  clear  from 
Fig.  3(a)  that  after  the  sudden  change  in  the  current  density, 
the  water  content  in  PEM  varies  significantly  during  the  initial 
period,  since  the  electro-osmotic  effects  decrease  the  water  con¬ 
tent  in  the  anode  and  increase  that  in  the  cathode.  After  about 
1  s,  the  water  content  near  the  cathode  side  approaches  to  the 
steady-state  value,  while  the  water  content  near  the  anode  side 
continuously  declines  with  time  until  it  approaches  the  steady 
state.  Comparison  of  corresponding  curves  between  Fig.  3(a) 
and  (b)  reveals  that,  as  humidification  parameter  k  is  increased, 
the  water  content  near  the  cathode  rises  and  quickly  approaches 
the  steady  state.  This  conclusion  can  be  plausibly  drawn  by 
noting  the  fact  that  a  larger  k  implies  a  higher  water  transport  per¬ 
formance  between  the  PEM  and  its  outside,  which  can  improve 
the  response  rate  of  water  transport. 

The  effects  of  humidification  on  both  cathode  and  anode  sides 
on  the  evolution  of  water  content  in  the  PEM  are  shown  in 
Fig.  4(a)  and  (b),  which  correspond  to  the  relative  humidity  of 
80%  and  60%,  respectively.  According  to  Fig.  4(a),  in  the  ini¬ 
tial  period,  the  water  content  near  the  anode  side  decreases  due 
to  the  electro-osmotic  effects,  while  the  water  content  increases 
near  the  cathode  side.  After  2  s,  however,  the  back  diffusion 
effects  begin  to  emerge  for  the  enlarged  water  content  difference 
between  the  cathode  and  anode  sides.  And  the  water  content  near 
the  anode  side  gradually  increases.  These  phenomena  are  more 
obviously  shown  in  Fig.  4(b),  for  the  relative  humidity  of  water 
vapor  on  the  cathode  and  anode  sides  equal  to  60%.  It  is  found 
that,  when  the  relative  humidity  on  the  anode  side  is  insufficient, 
the  water  content  distribution  in  PEM  varies  slowly  because  the 
back  diffusion  process  dominates,  thus  a  relatively  longer  time 
is  needed  for  the  system  with  a  lower  relative  humidity  to  reach 
the  steady  state. 

The  evolutions  of  the  water  content  in  the  PEM  at  relative 
humidity  of  water  vapor  of  80%  and  60%  on  the  cathode  side 
are,  respectively,  presented  in  Fig.  5(a)  and  (b),  with  the  relative 
humidity  of  water  vapor  on  the  anode  fixed  to  be  80%.  It  is 
clearly  seen  from  Fig.  5  that  when  the  water  content  near  the 
cathode  side  decreases,  the  transient  response  slows  down  near 


Fig.  3.  Effects  of  humidification  constant  on  the  evolutions  of  the  water  content 
in  PEM:  (a)  k  =  0.0 1  cm  s~ 1  and  (b)  k  =  0.00 1  cm  s“ 1 . 

the  cathode  side.  As  illustrated  in  Fig.  5(b),  when  the  relative 
humidity  of  water  vapor  on  the  cathode  side  decreases  to  60%, 
the  water  diffuses  to  the  cathode  side  in  the  initial  stage,  because 
the  water  content  near  the  anode  side  is  higher  relative  to  that 
near  the  cathode  side.  However,  as  the  operating  current  density 
increases  with  time,  the  water  generated  in  reaction  begins  to 
compensate  the  shortage  of  water  due  to  the  low  humidity  on 
the  cathode  side  in  the  initial  stage.  Then  the  back  diffusion 
effects  become  to  perform  for  the  large  water  content  difference 
between  the  cathode  and  anode  sides.  Consequently,  the  water 
content  near  the  anode  side  rises,  which  in  turn  causes  the  slow 
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Fig.  4.  Effects  of  humidification  on  both  cathode  and  anode  sides  on  the  evolu¬ 
tion  of  water  content  in  the  PEM:  (a)  80%  and  (b)  60%. 

response.  Therefore,  the  time  needed  for  the  water  content  in 
PEM  to  approach  the  steady  becomes  long. 

Fig.  6  presents  the  effects  of  the  thickness  of  PEM  on  the  evo¬ 
lution  of  water  content.  Due  to  the  electro-osmotic  effects,  with 
an  increasing  PEM  thickness,  the  water  content  near  the  anode 
side  of  the  PEM  continuously  rises,  while  the  water  content  near 
the  cathode  side  sustaining  drops.  As  their  content  difference 
grows  large  enough,  the  back  diffusion  effects  take  place.  If  the 
thickness  is  too  thick,  the  water  content  is  not  enough  for  the 
back  diffusion,  then  the  water  content  near  the  cathode  side  of 
the  PEM  declines.  If  the  thickness  is  small,  the  water  content  in 


Fig.  5.  Effects  of  humidification  on  the  cathode  side  on  the  transient  evolutions 
of  the  water  content:  (a)  80%  and  (b)  60%. 

PEM  varies  very  quickly  and  stably,  and  the  water  distribution 
becomes  more  uniform. 

Fig.  7(a)  and  (b)  presents  the  response  time  needed  for  the 
system  from  one  current  density  to  another  at  the  humidification 
constant  k  of  0.01  and  0.001  cms-1,  respectively,  with  the  rel¬ 
ative  humidity  of  water  vapor  being  100%  on  both  cathode  and 
anode  sides.  Additionally,  the  transverse  axis  denotes  the  initial 
current  density,  the  vertical  axis  denotes  the  response  time  for 
reaching  the  steady  state  from  one  current  density  to  another, 
and  the  three  curves  describes  three  different  switching  current 
densities.  The  comparison  of  Fig.  3(a)  and  (b)  shows  that,  when 
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Fig.  6.  Effects  of  thickness  of  PEM  on  the  evolution  of  water  content  in  PEM: 
(a)  Lm  =  0.015  cm  and  (b)  Lm  =0.005  cm. 

k  is  large,  the  response  time  for  switching  current  density  to 
1.0  A  cm-2  is  longer.  This  is  because  that  the  water  content  in 
anode  PEM  is  low  for  high  current  density.  When  k  is  large, 
however,  the  water  transport  between  PEM  and  its  outside  is 
enhanced,  so  that  the  distribution  changes  with  little  decrement, 
and  that  the  time  to  reach  the  steady  state  is  long. 

Fig.  8(a)  and  (b)  presents  the  response  time  for  reaching  the 
steady  state  from  one  current  density  to  another  with  the  relative 
humidity  of  water  vapor  on  the  cathode  and  anode  sides  being 
80%  and  60%,  respectively,  when  the  humidification  constant  k 
is  0.01  cm  s_1,  and  the  thickness  of  PEM  100  |jim.  It  is  found  in 


Fig.  7.  Effects  of  humidification  constant  on  the  response  time  required  for 
the  PEM  fuel  cell  reaching  the  steady  state  from  initial  current  density 
/o  =  0.1  A  cm-2  to  another  one:  (a)  k  =  0.01  cms-1  and  (b)  k  =  0.001  cms-1. 

Fig.  8(a)  that,  when  the  humidity  decreases,  the  response  time 
for  the  PEM  to  reach  the  steady  becomes  long.  Especially,  at 
a  low  operating  current  density,  the  electro-osmotic  effects  are 
weak,  and  consequently  the  water  content  in  PEM  is  mainly 
affected  by  the  back  diffusion  process.  Then  the  water  content 
in  anode  PEM  continuously  rises.  However,  at  a  high  operating 
current  density,  the  electro-osmotic  effects  are  enhanced,  and  the 
back  diffusion  effects  are  depressed,  so  that  the  response  time 
for  reaching  the  steady  becomes  short.  Additionally,  it  is  found 
in  Fig.  8(b)  that  the  response  time  needed  to  reach  the  steady 
is  longer  for  a  case  with  a  lower  humidity.  The  main  reason  is 
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Fig.  8.  Effects  of  humidification  of  both  anode  and  cathode  sides  on  the  response 
time  required  to  reach  the  steady  state  from  initial  current  density  /'q  =0.1  A  cm-2 
to  another  one:  (a)  80%  and  (b)  60%. 

that,  when  the  humidity  on  the  anode  side  is  very  low,  the  water 
content  difference  between  the  two  sides  of  PEM  is  so  large 
that  the  water  content  in  PEM  is  mainly  affected  by  the  back 
diffusion  process,  and  that  the  response  time  to  reach  the  steady 
is  enlarged.  As  a  conclusion,  at  a  low  humidity,  the  water  content 
in  PEM  is  insufficient,  and  the  response  time  for  the  fuel  cells 
to  reach  the  water  balance  increases. 

Fig.  9(a)  and  (b)  presents  the  response  time  to  reach  the 
steady  state  from  one  current  density  to  another  with  the  rela¬ 
tive  humidity  of  water  vapor  on  the  cathode  side  being  80%  and 
60%,  respectively,  when  k  is  0.01  cm  s_1,  the  thickness  of  PEM 


Fig.  9.  Effects  of  humidification  of  cathode  and  anode  on  response  time  required 
to  reach  the  steady  state  from  initial  current  density  Iq  =  0.1  A  cm-2  to  another 
one:  (a)  80%  and  (b)  60%. 

100  pun.  Comparing  with  Fig.  7(a),  the  response  time  to  reach  the 
balance  is  longer  with  switching  current  density  to  0.5  A  cm-2, 
while  the  response  time  with  higher  switching  current  density  is 
shorter  for  a  lower  initial  current  density.  It  is  mainly  due  to  the 
weaker  electro-osmotic  effects  at  lower  current  density.  When 
the  back  diffusion  effect  dominates  the  water  content  in  PEM, 
the  response  time  to  reach  the  balance  is  enlarged.  It  is  also 
found  in  Fig.  9(b)  that  the  response  time  for  the  water  content  in 
PEM  to  reach  the  steady  state  is  long  when  the  humidity  of  inlet 
water  vapor  is  insufficient.  If  the  humidity  of  inlet  water  vapor 
on  the  cathode  side  is  lower  than  that  on  the  anode  side,  the 
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response  time  to  reach  the  steady  state  at  higher  current  density 
is  shorter.  In  addition,  when  the  thickness  of  PEM  is  thickened, 
the  response  time  for  reaching  the  steady  state  increases  signif¬ 
icantly. 

4.  Conclusions 

The  present  work  investigates  a  one-dimensional  transient 
process  of  water  transport  in  PEM  fuel  cells.  The  diffusion  and 
electro-osmotic  effects  on  the  water  content  in  the  PEM  are 
explored,  with  only  the  diffusion  effects  are  considered  for  the 
vapor  liquid  transport  in  the  GDL  and  CL.  Corresponding  to  the 
results  of  analysis,  several  conclusions  are  drawn  as  follows: 

(1)  The  humidification  constant  k  plays  an  important  role  in 
the  water  transport.  When  k  is  large,  the  performance  of 
the  water  transport  from  PEM  to  its  outside  is  very  high, 
thus  the  water  content  in  PEM  is  uniform.  As  k  decreases, 
the  water  content  in  cathode  PEM  rises  due  to  the  electro- 
osmotic  effects,  and  the  excess  water  is  difficult  to  diffuse 
into  the  cathode  CL  and  GDL.  Similarly,  the  supplementary 
rate  of  the  water  content  in  anode  PEM  is  dependent  upon 
k,  thus  the  water  content  in  anode  PEM  is  lower  for  smaller 
k.  The  magnitude  of  k  also  affects  the  transient  behavior  of 
water  transport. 

(2)  When  the  humidity  conditions  on  the  two  sides  are  the  same, 
the  decrease  of  humidity  has  negative  effect  on  the  water 
content  distribution  in  PEM  and  the  response  time  to  reach 
the  steady. 

(3)  When  the  relative  humidity  on  the  cathode  decreases  to  80%, 
the  steady  distributions  of  water  content  in  PEM  under  var¬ 
ious  current  densities  have  little  difference  with  each  other. 
As  the  current  density  increases  from  0.1  to  1.0  A  cm-2, 
the  response  time  to  reach  the  steady  is  even  shorter  than 
that  with  100%  relative  humidity.  However,  when  the  rela¬ 
tive  humidity  in  the  cathode  decreases  to  60%,  the  steady 
water  content  in  PEM  remarkably  drops,  and  the  response 
time  to  reach  the  steady  state  becomes  significantly  longer. 
Therefore,  high  humidity  benefits  the  performance  of  fuel 
cells. 

(4)  When  the  PEM  is  thinned,  the  distribution  of  water  content 
is  more  uniform,  and  the  response  time  for  the  water  content 
to  reach  the  steady  is  apparently  shorter. 
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